Can a single gradientless light beam drag particles? 
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Usually a light beam pushes a particle when the photons act upon it. This is due to that the 
electric-dipole particle in the paraxial beam is considered. We investigate the scattering forces in 
non-paraxial gradientless beams and find that the forces can drag certain particles towards the 
beam source. The major criterion to be carried out to get the attractive force is the strong non- 
paraxiality of the light beam. The cone angle denoting the non-paraxiality has been investigated to 
unveil its importance on achieving dragging force. We hope the attractive forces will be very useful 
in nanoparticle manipulation. 



In spite of the well established theory, experiment and 
applications in particle-light interaction [1-5], there are 
many unknovifns yet to be explored in the particular field 
of attracting or separating molecules or nanoparticles by 
gradientless light beams. In the ordinary optical tweez- 
ers, the transfer of the particles along the 3D trajecto- 
ries can be achieved using the spatial light modulators 
[6]. In this case, the consecutive change of the computer- 
generated holograms for the tightly focused light beams 
produces the necessary effect, trapping of the particle due 
to the gradient of the field. 

It is well known that the light beam results in repelling 
force on the object. It is correct within electric-dipole ap- 
proximation and paraxial light beams, when the size of 
the particle is much less than the wavelength [2] . In this 
case, the force can still be manipulated to be a drag- 
ging force, provided that two oppositely directed beams 
are present to control the particle's position [7]. How- 
ever, such two-beam configuration (requiring exactly op- 
posite directions) is not necessary to achieve the attrac- 
tive force. As shown in Ref. [8] , two beams with different 
longitudinal wavenumbers can locally operate as a trac- 
tor beam. Another possibility to get attractive force is to 
use the gain media [9] , but the exotic media strongly con- 
fine the possible applications of the negative forces. In 
Ref. [10] it was suggested without proof that the nega- 
tive Poynting vector could be the reason of the attractive 
force. Though the negative energy flux density is not the 
reason of the particle drag towards the light source, this 
effect is quite possible in the non-paraxial beams. 

Small particles (hence, the dipole approximation) as 
a limiting factor are not necessary either. The particle 
can be large compared to the wavelength and positioned 
near the center of the beam. In this letter we use a 
non-paraxial gradientless beam (a vector Bessel beam) 
to induce the dragging force. The choice is owing to the 
following reasons: (i) long propagation distance, where 
beams weakly spread and keep their intensity profile; (ii) 
complex polarization structures of the beam; (iii) non- 
paraxiality of the beams denoted by the cone angle. We 
need such a general description of the beam, because it 



would be nice to know whether the beam should be ac- 
tually complex or we can use its simplified version. Field 
distribution in the form of the vector Bessel beam is not 
unusual. It can be generated as a mode of a circular 
fiber. However, the experimental generation of the non- 
paraxial beams is still challenging. 

The key idea is to demonstrate that a gradientless 
beam can exert a dragging force upon a particle. The 
beam being non-paraxial and gradientless is very impor- 
tant because the particle can thus be directly dragged 
toward the light source realizing the idea of the tractor 
beam. It should be noted that the Bessel beam recon- 
structs the field behind the particle [11]. This means that 
we can trap and drag several objects by a single beam. 

z-propagating electromagnetic Bessel beam in vacuum 
is described as [10] 

E = c'™^+'^^ f J™(gr)c2e, - ^ci(e, x b) + ^C2h) , 
V 1 1 J 

H = c'™^+'^- (^J™((?r)cie, + ^cib -t- ^C2(e, x b)^ ,(1) 

where fco = lo/c is the wavenumber in vacuum, q = 
fco sin a is the transverse wavenumber (a is the cone an- 
gle, see Fig. 1 (a)), (3 — \/k'^ — (f- is the longitudinal 
wavenumber, m is the beam order. 
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The typical intensity distribution |E|^ of the vector 
Bessel beam is illustrated in Fig. 1(b) and (c). The in- 
tensity is actually propagation invariant (z- independent) 
and azimuth invariant (</?- independent). It worth to draw 
attention to the coefficients ci and C2 of the vector Bessel 
beam (1). They provide the possibility to create the ap- 
propriately polarized propagation invariant field distri- 
bution and, in particular, can be associated with TE- 
polarized (c2 — 0) and TM-polarized (ci = 0) waves. 
In general, coefficients ci and C2 can be complex num- 
bers, what means that the TE and TM beams are phase 
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FIG. 1: (a) Light beam pushes a particle from position B' to 
position B while drags another particle from original position 
A' to new position A (it will keep dragging if the light keep 
shining). On the right, the gradientless bessel beam wavevec- 
tors lying on the cone are sketched, (b) and(c) denote the 
cross-sections of the intensity of such a beam propagating 
along z-axis (ci — 1, C2 = i, q/ko — 0.9, m = 1). 

shifted. The ordinary Bessel beams follow from Eq. (1) in 
paraxial approximation {q <C fco), but then we lose both 
possibility to set large qs and generate beams with com- 
plex polarization structure. Non-paraxial Bessel beams 
have already demonstrated a number of exciting proper- 
ties like local negative direction of the Poynting vector 
[10] and intensity transformation on reflection [12]. 

We assume that the spherical particle of the radius R is 
situated exactly at the center of the beam, but it is not 
a limitation, because the calculations can be made for 
the shifted particle as well. The incident vector Bessel 
beam is scattered by the particle according to the Mie 
theory. The scattered fields can be calculated, for exam- 
ple, using the matrix approach [13, 14], which is valid 
for any incident electromagnetic wave. So, the fields at 
the boundary of the sphere can be written as the sum 
of the incident and scattered fields: E = E™'^ -I- E'*'^ and 
H = H™'^ -I- H'*''. These fields are sufficient to compute 
the time-average electromagnetic force on the particle: 

/•TT /'27r 

F = / / {fnjR^smOdedifi, (2) 
Jo Jo 

where n is the outward normal, time-averaged Maxwell's 



stress tensor is 
f = -!-Re ( E eg) E* H H* - 1{\E\^ + IHp) ) . (3) 

OTT \ 2 / 

Here E E* defines the dyad, which in index form is 
(E (g) E*)ij — EiE*. We use the general formulae for 
calculations, but the scattered fields of the Rayleigh par- 
ticles can be written via the electric and magnetic dipole 
moments [3]. We are interested only in the component 
Fz of the force and look for the situations when < 0, 
i.e. when the particle is pulled by the beam. 

When the particle is located at the beam axis, angle (p 
for the Bessel beam (1) is simultaneously the azimuthal 
angle of the particle's spherical coordinates. Therefore, 
the dependence of the fields on (p keeps only as an ex- 
ponential exp{imip) and the incident and scattered fields 
contain the single term in the sum over the integer az- 
imuthal number. It is not the case for the particle shifted 
from the axis. Then we need to take into account also 
the terms with integer azimuthal numbers nonequal m. 

The structure of the beam strongly defines the possi- 
bility of the attractive force. However, the intensity of 
the light is not the criterium of the pulling. Even identi- 
cal intensity patterns (like in the insets of Fig. 2 (b) and 
(d)) result in the qualitatively different dependencies of 
the force. It is likely that the vector (polarization) struc- 
ture of the beam plays the leading role in the appearance 
of the attractive force. It should be noted that the nega- 
tive longitudinal component of the Poynting vector Sz is 
not responsible for the appearance of the negative force, 
because the large negative Sz exists for C2 — — i [10], 
while the force is positive in this case (Fig. 2(f)). 

Strong attractive force exists for non-phase-shifted su- 
perposition of TE and TM-polarized beams (c2 is positive 
real number as in Fig. 2(b) and (c)) and for 7r/2-shifted 
beams (c2 is imaginary number with Imc2 > as in Fig. 
2(e)). In the latter case, the negative force is broadband, 
low positive peaks being caused by the quadrupole terms 
(by the integer polar number ^ = 2 in the scattering se- 
ries). These two peaks at e = 8 and e = 10 exist in 
each subfigure and provide the positive force. The dipole 
term {I = 1) describes the features near e = 4. Due 
to the non-paraxial regime, and/or large particle size, 
and/or magnetic moments the force owing to the dipole 
term I — 1 can become attractive. 

Thus, the Bessel beams with C2 = 1 and C2 = i are 
preferred, but they are quite different. In Fig. 3 we show 
the force acting on the particle for C2 = 1 (subfigures (a) 
and (c)) and C2 — i (subfigures (b) and (d)). Asymmetry 
in the case of C2 = 1 implies that, in order to induce the 
attractive force the permittivity should be greater than 
the permeability, i.e. the excited electric moments play 
the main part. In the symmetric case (Fig. 3(b) and 
(d)), both electric and magnetic inputs are of the same 
order and the force has a minimum at the impedance 
matching line e — fi. It is interesting that the region of 
the negative force is really great. The value of the force 
can be enhanced by adjusting the radius of the sphere 
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Dielectric permittivity e 



FIG. 2: Force Fz on the spherical particle {koR — 1, n — 3) versus dielectric permittivity e for (a) C2 — 0, (b) C2 — 1, (c) 
C2 = 2, (d) C2 = —1, (e) C2 = i, and (f) C2 — — i. In the insets, the field intensities are demonstrated (grey circle stands for the 
particle). Beam parameters are ci — 1, q/ko = 0.9, and m — 1. 



(see Fig. 4(a)). From the figure Fig. 3(d) we can guess 
that the force can be attractive even for = I (it is out 
of the range of the plot). It is actually the case as can 
be seen from Fig. 4(a). So, the dragging force can act 
even on the common glass particles, if it has an appro- 
priate size. The maximizing of the absolute value of the 
negative force using the impedance matching condition 
substantially broadens the range of the negative force. 
The only strong limitation on the vector Bessel beams is 
their non-paraxiality. As shown in Fig. 4 (b), the varia- 
tion of the parameters does not substantially reduce the 
cone angle a — arcsin(g/fco) (for q ~ 0.9 the angle a 
equals ss 64°). The beam should have small longitudinal 
wavenumber /3 to be able to create the attractive force. 

From Fig. 4(a) we can suggest that even for the 
Rayleigh magnetodielectric particles the force can be less 
than zero. In fact, the mechanism behind the attractive 
force can be revealed from the theory of magnetoelectric 
Rayleigh particles [3] . Using the fact of the small imag- 
inary part of the polarizability in comparison with the 
real part, we derive 

(F.) = ^ (lm(ae)|Ep +ImK0|H|2) 

-^Re(ae)Re(a„)Re(P,), (4) 

where Pz — ez(E x H*), ae and a„i are the electric 
and magnetic polarizabilities of the particle, respectively. 
From Eq. (4) it is evident that the negative force is 
feasible due to (i) small longitudinal wavenumber /3 of the 
light beam, (ii) existence of the magnetic dipole moment. 
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FIG. 3: Diagram of the force Fz as function of the dielectric 
permittivity e and magnetic permeability ^ for (a), (c) C2 = 1 
and (b),(d) C2 = i. In the bottom subfigures, only negative 
values of the force is shown. Parameters: koR — 1, ci = 1, 
q/ko — 0.9, m = 1. 

and (iii) large positive Poynting vector (quantity Re(Pz)). 

If /3 is not small, what is the case of the small spherical 
particles {ae = a^P^ + i2/cQQ;e°''^/3, where a^P^ = a^(e — 
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Cone angle a, deg. 

FIG. 4: Dependence of the z-component of the force on 
(a) the radius of the spherical particle fcoi? (g/fco = 0.9) 
and (b) cone angle a = arcsin(5/fco) (ci ~ i). A multime- 
dia file is created to show the significance of the cone angle 
(non-paraxiality effect) and particle's radius (size effect), in 
order to achieve dragging forces (see Media 1. Hyperlink: 
http://www.ece.nus.edu.sg/stfpage/eleqc/force.gif). Param- 
eters: ci = f, m = f. 



l)/(e -l- 2), and the similar expressions for the magnetic 



polarizability) in the field of the plane wave (/? = 1 and 
|E|2 = |H|2 = Re(P^) = 1) we write 

fc|a6 f je-lf (£-1)(a^-1) \ 
^ 3 V(^ + 2)2"^(/i + 2)2 (£ + 2)(^ + 2);- 

The force cannot be attractive for any dielectric per- 
mittivity and magnetic permeability under the illumina- 
tion of a plane wave. 

For non-magnetic small particles (am — 0), the force 
is always pushing (it follows straightforwardly from Eq. 
(4)). Nevertheless, non-magnetic non-Rayleigh particles 
can be attracted owing to the contributions of the higher- 
order electrical moments. 

Eq. (4) claims the Poynting vector should have large 
positive value. From [10] we notice that the largest value 
of the Poynting vector is achieved for C2 = i and this 
fact is in full agreement with the calculations in Fig. 2. 
However, we should be careful in generalization of the 
results for the Rayleigh particle to the non-Rayleigh ones. 
Indeed, the Poynting vector distributions are the same for 
C2 — 1 and C2 — —1, but the forces dramatically different 
in the case kgR = 1 as it is observed in Fig. 2 (b) and 
(d). So, Eq. (4) does not contradict the existence of the 
attractive force for magnetodielectric Rayleigh particles. 

In conclusion, we have demonstrated the possibility 
of the attractive forces from a gradientless light beam. 
Compared to the previously reported tractor beams [8], 
we have the broadband negative force and do not need 
two opposite beams with different wavenumbers [7] or 
one gradient beam [6]. The proposed gradientless beam 
with the proper manipulation from the non-paraxiality 
can lead to a lasting dragging force exerted on particles 
which can be large or small and made of any materials. 
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